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We present first results, for heavy to moderate quark masses, of a study of thermodynamic properties of 1-flavor 
QCD, using the multiboson algorithm. 



> 

m 
O 

On 
O 
\D 



X 



1. Introduction 

Our understanding of QCD thermodynamics is 
still very incomplete. Simulations of the 0-flavor 
(quenched) lattice theory are just starting to pro- 
duce results in line with continuum expectations. 
The situation in the 2- and 3-flavor cases is still 
rather murky (for an excellent review, see or 
this volume). Additional information for different 
numbers nt of quark flavors is clearly desirable to 
build a coherent picture such as proposed in Q 
near the chiral limit. The study of QCD with ex- 
otic ft/ has focused on many flavors Perhaps 
because of algorithmic difficulties, 1-flavor QCD 
has been largely ignored, although it has partic- 
ularly interesting properties. 

It is impossible to make a tt — a doublet with 
only 1 quark species, so there is no chiral symme- 
try to speak of at vanishing quark mass. The 
low-lying meson spectrum consists of the "/o" 
(qq) and the "77"' (975(7). The 77' correlator 
G n >(x,y) = \D-\x,y)\ 2 - D- 1 {x,x)D- x {y,y) 
contains 2 pieces, only the first of which, anal- 
ogous to the pion correlator, becomes massless 
(constant) at m q = 0. Thus the chiral phase tran- 
sition and the associated Goldstone boson are ab- 
sent in the 1-flavor theory . In the heavy quark 
regime, on the other hand, the effect of dynami- 
cal quarks is still analogous to that of an exter- 
nal magnetic field on the equivalent Potts system. 



Just like in the 2-flavor case, this external field 
weakens the phase transition, possibly turning it 
into a crossover. Thus, as the quark mass is de- 
creased, one might expect the first-order decon- 
fincment transition to become second-order, with 
an end-point at some finite quark mass. For any 
lighter quark, only a crossover would be seen. 

Of course, testing this scenario requires care- 
ful finite-size scaling studies. As a first step, we 
present here preliminary results obtained on an 
8 3 x 4 lattice. We use Wilson fermions, and try 
to follow the critical line in the (k, (3) plane, start- 
ing from k = 0. 

2. Algorithm 

The standard algorithm to study any number 
of flavors is the R-algorithm Q], where the ef- 
fect of the fermionic determinant det_D is replaced 
by discretized Brownian noise. This discretiza- 
tion entails step-size errors, so that a delicate 
and costly extrapolation to zero step-size must 
be performed. Instead, we use here the multibo- 
son method In its non-hcrmitian variant, it 
can simulate 1-flavor QCD 0. 

Approximating D~ x by a Chebyshev polyno- 
mial P 2n {D) = c 2n ]Xl =1 {D - z k )(D - Zk), and 
noticing that det(Z? — Zk) = det75(Z? — 2^)75 = 
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det(Z? — Zk)', one obtains 
det£> re det~ 1 P 2n {D) = \detT n (D)\' 



(1) 
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with 



T n {D) 



Zk) 
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The right-hand side of ([[]) can then be replaced 
by a Gaussian integral over n bosonic fields 4>k- 
The action is S to t = S gaugc + J2k=i \(D - z k )(/>k\ 2 ■ 
Because (Q) is an approximation, a corrective 
Metropolis test is added, just like in the 2- 
flavor algorithm |J. At the end of each "tra- 
jectory" {U, (j)} — > {U',(j)'}, the configuration is 
accepted with probability min(l, q-v (W w-i)r]^ 
where W = PU^T" 1 {D')\[T m {D)T-\D))-\ 
This algorithm samples the dis- 
tribution e _s,sauBO det _1 P2m(-D)- With negligible 
overhead, m can be taken very large (m > 3n 
in our case), so that the measure is effectively 
e - s «a^det -1 P 0O (D). 

Now the domain of convergence S of the Chcby- 
shev approximation P2m(z) ~ 1/z in the com- 
plex plane is bounded by an ellipse dS centered 
at (1, 0) which goes through the origin. If z € S, 
zPim(z) = l+(D(e~ am ): the approximation is ex- 
ponentially convergent. Otherwise, |zP2 m (,z)| — > 
oo as ttl — > oo . This means that configurations 
for which eigenvalues of D fall outside S will al- 
ways be rejected by our algorithm. The measure 
we sample is e gaugo detZ? 8(spectrum(D) £ S), 
where O(u) — 1 or if u is true or false. 

We take for dS a circle, which contains the 
spectrum of D = 1 — kM at (3 — oo for any 
k < k c ~ 1/8, and at (3 — for any k < k c — 1/4 
in the thermodynamic limit. Thus we expect it 
to contain the spectrum of D also at interme- 
diate [3 for k < k c (/3), in the thermodynamic 
limit. Finite- volume fluctuations cause some "ex- 
ceptional configurations", with eigenvalues out- 
side S, to be included in the exact partition func- 
tion. Among those are configurations with nega- 
tive real eigenvalues, for which det-D < 0: they 
create the notorious "sign problem" associated 
with non-degenerate Wilson quarks. Our algo- 
rithm discards such configurations 

Therefore, even though it recovers the correct 
measure in the thermodynamic limit, our algo- 



rithm introduces a bias on a finite lattice. This 
bias however can be viewed as a finite-size effect: 
it only appears if the inverse quark mass becomes 
too large compared with the lattice size. For the 
relatively heavy quarks we consider here, this bias 
is completely negligible. 

The implementation of our algorithm on the 
Cray-T3D is very similar to the 2-flavor case ||. 
Here too, even-odd preconditioning can be imple- 
mented through a simple redefinition of the zeroes 
Zk ■ Half as many fields <pk are needed for the same 
quality of the approximation (0). The Metropo- 
lis acceptance is shown in Fig.l, to be compared 
with the 2-flavor case. 



2 Note that the R-algorithm would include them in the 
partition function, but with the wrong sign. 



Figure 1. Acceptance versus number of auxiliary 
fields for 4 4 , (3 = 0, k = 0.215. Circles are n f = 1 
data. Diamonds are rij = 2 data. 

3. Results 

For a fixed k, we tried to scan (3 and iden- 
tify "/3 c (/v)". On our rather small lattice, differ- 
ent criteria for criticality will give slightly differ- 
ent effective (3 C 's. We tried to choose a criterion 
which reduced finite-size effects, at least in the 
quenched case. Susceptibility measurements, of 
the Polyakov loop especially, had large statisti- 
cal errors and were difficult to exploit. The most 
useful observables for us were: 
i) the histogram of the norm of the Polyakov loop, 
which shows 2 peaks near criticality. This is safe 
only if many tunneling events occur between the 
2 phases. Fortunately that is the case with our 
algorithm, as shown in Fig. 2 for n = 0.10. The 
work per sweep is roughly equivalent to 6 molec- 
ular dynamics steps of the R-algorithm, which 
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shows that the multiboson approach is rather ef- 
ficient to decorrelate the Polyakov loop, a global 
observable. 

ii) the deconfinement ratio p = 3/2p— 1/2, where 
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Figure 2. History of the norm of the Polyakov 
loop for f3 = 5.67, k = 0.1, L = 8 3 x 4. The 
work per trajectory is equivalent to ~ 6 molecular 
dynamics steps. 

p is the probability for the complex trace of the 
Polyakov loop to fall within ±20° of a Z3 axis, p 
goes from to 1 as the Z3 symmetry is broken. A 
preliminary quenched study showed us that iden- 
tifying (3 C as that for which p = 75 — 80% gave 
very small finite-size corrections. 

We then used these 2 observables to deter- 
mine (3 C for k = 0.05,0.10,0.12,0.14. The num- 
ber of auxiliary fields cj>k was 8, 16, 24, 32 respec- 
tively. Representative histograms of the norm of 
the Polyakov loop are shown in Fig. 3, to compare 
with the quenched case. A 2-peak structure re- 
mains visible up to the highest k. One can also 
clearly see the breaking of the Z3 symmetry by 
the dynamical quarks: the peak corresponding to 
the confined (disordered) phase moves to the right 
as k increases. This is evidence for the expected 
weakening of the phase transition, which should 
be confirmed by simulations on larger spatial vol- 
umes and by finite-size scaling. 

Based on Fig. 3 and on the measurement of the 
deconfinement ratio, our estimates for [3 c (k) are 
given in Table 1, together with published results 
for nf = and 2. One sees little departure so far 
from a linear behavior in nj: the ordering effect 
of 1 flavor of quarks is about 1 /2 that of 2 flavors. 
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5.69254(24) 


0.05 


5.685(5) 




0.10 


5.67(1) 




0.12 


5.630(15) 


5.58(2) 


0.14 


5.59(2) 


5.46(2) 



Table 1: /3 c (k) for the quenched ]8j, one- and two- 
flavor theories 
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Figure 3. Histograms of the norm of the Polyakov 
loop for parameters near C . 



